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The ISM in spiral galaxies: can cooling in spiral shocks 
produce molecular clouds? 



00 

o 

O 



(N 

Oh 
6 



> 

(N 
m 

O 

oo 
o 



X 



C. L. Dobbs, 1 * S. C. 0. Glover, 2 P. C. Clark 3 & R. S. Klessen 3 

1 School of Physics, University of Exeter, Stocker Road, Exeter, EX4 4QL 

2 Astrophysikalisches Institut Potsdam, An der Sternwarte 16, D-144S2, Potsdam, Germany 

3 Institut fur Theoretische Astrophysik, Universitdt Heidelberg, Albert- Ueberle-Str. 2, Heidelberg, Germany 



26 June 2008 



ABSTRACT 

We investigate the thermodynamics of the interstellar medium (ISM) and the forma- 
tion of molecular hydrogen through numerical simulations of spiral galaxies. The model 
follows the chemical, thermal and dynamical response of the disc to an external spiral 
potential. Self-gravity and magnetic fields are not included. The calculations demon- 
strate that gas can cool rapidly when subject to a spiral shock. Molecular clouds in the 
spiral arms arise through a combination of compression of the ISM by the spiral shock 
and orbit crowding. These results highlight that local self-gravity is not required to 
form molecular clouds. Self-shielding provides a sharp transition density, below which 
gas is essentially atomic, and above which the molecular gas fraction is > 0.001. The 
timescale for gas to move between these regimes is very rapid (^1 Myr). From this 
stage, the majority of gas generally takes between 10 to 20 Myr to obtain high H2 
fractions (> 50 percent). These are however strict upper limits to the H2 formation 
timescale, since our calculations are unable to resolve turbulent motions on scales 
smaller than the spiral arm, and do not include self-gravity True cloud formation 
timescales are therefore expected to be even shorter. 

The mass budget of the disc is dominated by cold gas residing in the spiral arms. 
Between 50 and 75 percent of this gas is in the atomic phase. When this gas leaves the 
spiral arm and drops below the self-shielding limit it is heated by the galactic radiation 
field. Consequently, most of the volume in the interarm regions is filled with warm 
atomic gas. However, some cold spurs and clumps can survive in interarm regions for 
periods comparable to the interarm passage timescale. Altogether between 7 and 40 
percent of the gas in our disc is molecular, depending on the surface density of the 
calculation, with approximately 20% molecular for a surface density comparable to 
the solar neighbourhood. 
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1 INTRODUCTION 

Most stars form in giant molecular clouds (GMCs), and 
so understanding how GMCs form and disperse is a ma- 
jor goal of star formation research in a Galactic con- 
text. GMC formation is likely to be strongly dependent 
on the nature of the interstellar medium (ISM). Two 
prominent theories that have previously been advanced are 
that GMCs form by large scale magnetic and/or gravita- 
tional jiistabJUtiesJn t hermally or turbulentl y supported 
gas jElmcgrccn] 1 19791: fealbus fc Cowiel [l985: Elm egreenl 
11994 IChou et alj|200d : iKim fc Ostriken 120061 "). or by the 



collisional build up of smaller clouds or clumps (e.g 



Field fc Saslaw 19651: Levinson fc Robert si 1981 : Tomisakal 



19841 : iRoberts fc Stewartlll987l : iKwan fc Valdeslll987l ). The 

difference between these two scenarios is that in the first 
the ISM behaves essentially as a fluid, whilst in the latter 
the behaviour is more ballistic. The first case resembles the 
warm, low density component of the ISM, whilst the second 
appears more consistent with the cold HI or even molecular 
component which is predominantly situated in clumps. In 
the majority of galaxies, the ISM may be expected to con- 
sist of both of these components, with both i nstabilities and 
collis i ons contribu ting to GMC formation l|Zhang fc Sond 
ll999l : lDobbsll2008l ). 
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A major limitation of many galactic simulations is 
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that the ISM is assumed to be isothermal, either warm 
dKim fc Ostrikerl 12002): ISlvz et all 120031; IChakrabarti et all 



120031: iLi et alJ 120051. 120061: JShettv fc Ostrikerl l2006h cold 
HI (Dobbs & Bonncll 2006), or a two phase medium 
l|Dobbs fe Bonnelll l2007fl . In this paper we include a com- 
prehensive treatment of the thermodynamics of the ISM in 
global simulations, to investigate whether cooling in spiral 
shocks produces a clumpy, molecular phase even in the ab- 
sence of self gravity. The only other calculations of heating 
and cooling in grand design galaxies are presented in the 
recent paper bv lWadal (|2008l ). who uses an energy equation. 
Here it is possible to see how the nature of the shock changes 
with time as the gas changes from a warm smooth distribu- 
tion to cold clumps. A similar approach has been adopted for 
simulations without a spiral potent ial, which show complex 
thermal and den s ity distributions llWada fc Normanll200ll : 
IWada et al.ll2002l ; pTasker fc Brvanll2006l ). 

Another key issue related to GMC formation is the de- 
termination of the timescale over which these clouds are 
formed. Although there is no generally accepted definition 
of the formation timescale, one plausible measure is the time 
required for the gas to become fully molecular. To address 
this, simulations ideally need to include the chemistry of H2. 
The formation of H2 ha s recently been modell ed in simula- 
tions of dwarf galaxies l|Pelupessv et al. 2006), but only a 
very simplistic treatment has been ap plied to spiral galaxies 
dWada et al.ll2000l : lDobbs et alj|2006h . 

Simulations which include a treatment of the thermo- 
dynamics of the ISM typically focus on much smaller length 
scales than th ose we study here (see , for instance, the re- 
cent review bv lHennebelle et al.ll2007h . Simulations of shock 
compressed layers show the formation of structure in the 
cold neutral medi um (CNM), and the g eneration of low lev- 
elsof turbulence |Kovama & Inutsuka 2002; Hcitsch~et al.l 



2005; IVazquez-Semadeni et al.l 120061 : iHennebelle fc Auditl 
20071 ). Even without self g ravity, densities in excess 



of 1000 cm 6 are produced (|Audit fc Hennebelld [2005). 
The inclusion of self gravity leads to star formation 
in the density enhancements from the initial turbulence 
(|Vazquez- Scmadcni et al. |l2007l ). In addition to the thermo- 
dynamics, I Glover fc Mac Low! (|2007bl ) also follow the evo- 
lution of molecular hydrogen, They find that in a turbulent 
environment, H2 formation takes place over timescales of a 
few Myr, suggesting that molecular clouds may form very 

quickly. 

In previous work l|Dobbs et al.l l2006), we showed that 
molecular clouds c an form from cold gas by agglomeration 
in a spiral shock l|Dobbsl [20081 ). However these isothermal 
calculations assumed that either cold gas is prevalent in the 
ISM, or that warm gas cools very quickly in spiral shocks. 
Here we investigate whether those calculations are plausible, 
by including cooling and heati ng of the ISM. We in c lude th e 
thermodynamics described in iGlover fc Mac Low! (|2007al ). 
which takes into account the main processes responsible for 
cooling and heating of the ISM, rather than assuming a 
prescribed relationship between density and temperature. 
There are essentially two possible outcomes. If cooling is 
very rapid in the spiral shock, cold clumps are likely to form 
and the agglomeration of these clumps into GMCs is then 
a possibility. In this scenario, we are also interested to see 
whether the gas remains cold in the interarm regions, and 
whether spurs can be produced. Alternatively, cooling may 



not be particularly efficient, and little structure is produced 
in the spiral arms. In this case, another mechanism would 
be required to form GMCs. 



2 CALCULATIONS 

We perform these calculations using smoothed particle hy- 
drodynamics (SPH), a Lagrangian fluids code. The code is 
based on an original version by Benz (|Benz et al.l Il990l ), 
but has since been subject to significant modifications. The 
most substantial changes include the a ddition of sink par- 
ticles <lBatelll995l) a nd magnetic fields l|Price fc Monaghanl 
120051 : iPrice fc B ate 2007), although they are not used in the 
current paper. The code also includes individual timesteps 
and variable smoothing lengths. The density and smoothing 
lengths are solved iteratively according to 
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(1) 



l|Price fc MonagharJ|2007h . where h is the smoothing length, 
p the density, m the mass of the particle, and v is a di- 
mensionless parameter set to 1.2 in order that each particle 
has ~60 neighbours. Artificial viscosity is included to treat 
shocks, using the standard pa rameters, a = 1 and /3 = 2 
|Monaghan fc Lattanzidll985l '). 

2.1 Galactic potential and initial distribution 

We model a gaseous disc, and represent the stellar com- 
ponent of the galaxy by an external potential. The exter- 
nal potential includes a stellar disc, halo and incorporates 
a 4 arme d spiral compone n t. Th e potential is described 
in full i n iDobbs fc Bonnelll J2 OO6). The spiral component 
is from ICox fc Gomezl (120021 ) and has a pattern speed of 
2 x 10~ 8 rad yr" 1 , and a pitch angle of 15°. The amplitude 
of the stellar spiral perturbation is 1.1 x 10 12 cm 2 s -2 . 

Particles are initially distributed randomly in a torus 
with radii 5 kpc < r < 10 kpc. The gas particles are as- 
signed circular velocities according to the disc potential, and 
in addition, a Gaussian velocity dispersion of 6 km s _1 is im- 
posed (which constitutes the z component of the velocities) . 
The spiral perturbation emerges with time as the simulation 
progresses. Initially all the gas has a temperature of 7000 K 
and a scale height of 400 pc. We perform 3 simulations, with 
surface densities of 4, 10 and 20 Mq pc -2 . Although most 
galaxies do not have average surface densities as low as 4 Mq 
pc -2 , the lowest density calculation would correspond to the 
outer regions of galaxies. The total surface density in our 
calculations also includes helium, which has an abundance 
(by number) one-tenth that of hydrogen, i.e. n(He) = O.ln 
where n(He) is the number density of helium atoms and n 
is the number density of hydrogen nuclei (in all forms). 

We neglect both self gravity and magnetic fields in these 
calculations. These are discussed in p revious calculations us- 
ing an isothermal equation of state dKim fc Ostrikerl [20021 : 
IShettv fc Ostrikerll2006l : IDobbs fc Pricell2008l : |Pobbsll2008l ). 
Rather, we are interested in whether spiral shocks alone are 
sufficient to cool gas to temperatures of sC 100 K and pro- 
duce significant amounts of molecular gas. In the next sec- 
tions we describe the H2 chemistry in these simulations, and 
outline the main thermodynamic processes included. 
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2.2 Formation of H 2 



The abundance of mol ecular gas is c ompu ted according to 
the expression given in iBergin et all (|2004h . 



dn(U 2 
dT 



= R gr (T)nn(H) - [C CT + CdJaa(iV(H 2 ), A v )]T»(H a ), 

(2) 

where n(H) and n(H2) denote the number densities of H 
and H2, respectively, n = n(H) + 2n(H2) is the total num- 
ber density of hydrogen nuclei, AT(H 2 ) is the column den- 
sity of molecular hydrogen, Av is the visual extinction and 
T is the temperature. The formation rate of H2 on grains 
is R gr (T)nn(B), where R gr = 2.2 x 10~ 18 ST ' 5 cm 3 s _1 . 
The formation rate assumes that the num ber density of dust 
grain s is proportional to the total density (|Hollenbach et ahl 
Il97ll ). with the ratio of the dust grain number density to 
the hydrogen nuclei number density n being incorporated in 
Rgr- The parameter S is the efficiency of H2 formation on 
grains and is assumed to be a constant value of 0.3. In prac- 
tice, S depen ds on both the gas temperature and the grain 
temp erature (jHollenbach fc McKeelll989l :IC azaux fc Tielend 
|200J). At high temperatures, this is unimportant, as the 
density of the warm gas in our simulations is in any case too 
low to prevent immediate dissociation of alm ost all of the 
H2 formed on the grains |Dobbs et al.ll2006h . In cold gas, 
on the other hand, we may underestimate the H2 formation 
rate somewhat. 

The H2 is primarily destroyed by photodissociation. 
The photodissociation rate in these calculations is given by 

U SS (N(R 2 ),A V ) = / shic id(Af(^2))/dustCdiss(0) (3) 

( Drai ne fc Bertoldj[l996h . where £diss(0) is the photodissoci- 
ation rate for unshielded H2, /shield is a factor accounting for 
the effects of H2 self-shielding, and /dust is a factor account- 
ing for the effects of dust absorption. For the unshielded 
photodissociation rate, we take Cdiss(O) = 4.17 x 10 -11 s _1 , 
as appropriate for an ultraviolet radiation field equivalent 
to the galactic average from B O stars. For /shield, we use the 
following fitting function from iDraine fc Bertoldil (|1996| ) 

, ,. rm u 0.965 0.035 

./shield (,JV(H 2 J ) — 7— + 7— wTT , s 

(l + x/3) 2 (l + :r) u - 6 (4) 
x exp[-8.5 x 10~ 4 (1 + a;) - 5 ], 

where x = 7V(H 2 )/5x 10 14 cm -2 , and where we have assumed 
a typical Doppler br oadening parameter b = 3 km s _1 . For 
/dust, we again follow [Draine fc Bertoldil l| 19961 ) and assume 
that 



/dust — e 



(5) 



where t^iooo is the optical depth due to dust at a wavelength 
A = 1000 A, given by r dilo0 o = 2 x lO" 21 ^ + 2Ah 2 ). For 
gas with Ry = Av/E(B — V) = 3.1, as is typical of the 
diffuse ISM, r^iooo and the visual extinction Av are related 
by T dj iooo = 3.74Ay. 

In addition to photodissociation, H2 is subject to de- 
struction by cosmic rays. We assume a constant cosmic ray 
ionisation rate of £ cr = 6 x 10 -18 s _1 , although this term 
is minimal compared to the photodissocation rate, unless 
iV(H2) or Ay are very large. 

During the simulation, the fraction of molecular gas is 
always non-zero; the H2 is never completely photodissoci- 
ated. When the H2 fraction is very small, the photodis- 



sociation rate is large, and the requirement that n(H2) + 
Afchem x dn{R2)/dt > forces us to use a chemical timestep 
Atchem that is small compared to the other timescales in the 
problem. To accommodate this without compromising the 
overall efficiency of the code, we subcycle the evolution of 
the H2 abundance, splitting each hydrodynamical timestep 
into as many chemical su bsteps as are required (see also 
Gl over fc Mac Lowll2007a! ). The abundance of H 2 is used to 
determine the cooling rate, and so is consistently coupled 
with the evolution of the energy during the simulation. 



2.2.1 Calculating the H and H2 column densities 

The H2 photodissociation rate is strongly coupled to both 
the H2 column density (through /shield) and the total col- 
umn density of hydrogen nuclei, JVtot = ^V(H) + 2iV(H 2 ) 
(through /dust). To calculate these column densities, we use 
the simple estimate that the column density (of H2 or hy- 
drogen nuclei, as appropriate) is just the local density times 
the typical distance to a B0 star, a constant length scale, Z p h- 
The average distance to a B star is likely to be between 20 
and 50 pc depending on which spectral types are include d 
|Mihalas fc Binnevlll98ll ; lReedll200(]| ; lMafz-Apellanrzll200ll ). 
We therefore test 3 different length scales in these simu- 
lations. Our fiducial value is l p h = 35 pc, but we perform 
calculations with 15 and 100 pc, to estimate the minimum 
and maximum amo unts of H2 that may be produce d. In pre- 
vious calculations (|Dobbs et al. 2006; Dobbs 2008) we took 
lph = 100 pc. 



2.3 Thermodynamics 

To follow the thermal evolution of the gas in our simulations, 
we use a model for the h eating and cooling o f atom ic and 
molecular gas taken from Gl over fc Mac Low! (|2007aD . This 
model incorporates cooling from the fine structure emission 
lines of C + , O, and Si + , ro- vibrational emission from H2, 
gas-grain energy transfer, and recombination on grain sur- 
faces. In hot gas, additional cooling comes from the colli- 
sional dissociation of H2, the collisional ionisation of atomic 
hydrogen, and emission from atomic resonance lines (e.g. 
Lyman-a) and bremsstrahlung. Heating is provided primar- 
ily by photoelectric emission from dust grains, with addi- 
tional contributions coming from cosmic ray ionisation, H2 
photodissociation, and the pumping of excited vibrational 
states of H 2 by UV excitation or during the formation pro- 
cess of the molecules. Full details of th e rates adopted for 
these v arious processes can be found in lGlover fc Mac Low! 
(|2007ah . 

In dense gas, dust extinction sign ificantly redu c es the 
photoelectric heating rate. We follow IBergin et al.l (|2004l ) 
and assume that the heating rate scales with the visual ex- 
tinction Av as 



T pe (Av) = r pc (0)e 



-2.5A V 



(6) 



where r pc (0) is the photoelectric heating rate in optically 
thin gas, and where Av is calculated as in section [2~2l above. 

Aside from our simplified treatment of the dust extinc- 
tion, the main approximation that we have made is our 
neglect of the effects of emission from atomic carbon and 
silicon (which are not tracked in our chemical model) and 



4 C. L. Dobbs, S. C. Glover, P. C. Clark, R. S. Klessen 



Model 


E 

(M e pc- 2 ) 


'ph 
(PC) 


No. particles 


Mp ar t 

(Mq) 


1 


4 


35 


8 xlO 6 


125 


2 


10 


35 


8 xlO 6 


312 


3 


20 


35 


8 xlO 6 


625 


4 


10 


15 


8 xlO 6 


312 


5 


10 


100 


8 xlO 6 


312 


6 


10 


35 


1 xlO 6 


2500 


7 


10 


35 


4 xlO 6 


625 



Table 1. Table listing, for each simulation, the surface density, 
photodissociation length scale (/ p h), the number of particles and 
the particle mass (M part ). 

from CO. This simplification means that the temperature 
of the dense, highly molecular gas in our simulations will 
not be computed entirely correctly. However, at the densi- 
ties probed by our current simulations, the errors involved 
should be small, since C + cooling is comparable in effective- 
ness to CO cooling, as the significant quantities of cold gas 
produced in our simul ations readily demonstr ate (see e.g. 
section 13.1.11 see also [Gi over fe Jappsenll2007l ). Moreover, 
we would expect the dynamics of this cold, dense compo- 
nent to be dominated by turbulent motions, rendering its 
precise temperature of limited importance. 



3 RESULTS 

Details of the different calculations performed are shown in 
Table 1. We adopt 3 different surface densities, since the 
degree of collisional cooling in the disc is expected to de- 
pend predominantly on the density. As described in Sec- 
tion 2.3.1, three different length scales are used to calculate 
the photodissociation rate, which determines the molecular 
gas abundance. The next sections describe the structure of 
the disc for the different surface densities, and the resulting 
thermal distributions. The evolution of H2 is discussed in 
Section 3.2, including the amount of gas in the disc which 
becomes molecular and the timescales for H2 formation. For 
the simulations in the main part of the paper, we use 8 
million particles, but we also performed simulations with 1 
and 4 million particles to investigate the dependence of our 
results on resolution (see Appendix). 

3.1 Structure of the disc 

The galactic disc is shown for the different surface densi- 
ties in Figs 1 and 2. In all cases, the disc exhibits a large 
degree of substructure. The morphology of the disc is sim- 
ilar to that seen in previous isothermal calculations of cold 
gas IjDobbs fc Bonnellll2006l ). This is because, as will be dis- 
cussed in Section 3.1.1, a substantial amount of cold gas is 
formed in all of these calculations, and dense features such 
as clumps and spurs are produced in the cold gas. 

Figure 1 shows the structure of the disc for our standard 
model, with E = 10 Mq pc~ 2 and Z p h = 35 pc, at 3 times 
during the evolution of the disc. At the earlest time, there is 
not substantial substructure, but small spurs are beginning 
to emerge from the spiral arm (the transition between gas 
that has passed through a strong shock and gas which has 
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Figure 1. The evolution of the galactic disc is shown for the 
standard case, where E = 10 Mq pc -2 and £ p h = 35 pc. The 
disc contains a substantial degree of su bstructure, in agreemen t 
with previous isothermal calculations llDobbs fc Bonncll 20061) . 
Gas cools and becomes dense in the spiral arms, leading to the 
agglomeration of clumps into larger structures which shear off the 
spiral arms to become spurs (lower panels). 

not is evident: up to a radius of between 7 and 8 kpc, all 
the material has passed through at least one shock) . At this 
stage, there is cold gas in the spiral arms, but little sub- 
structure has yet developed. By 240 Myr, the structure in 
the disc is very similar to previous isothermal calculations 
ijDobbs fc Bonnelllliool ). with spurs extending into the in- 
terarm regions. The spiral shock causes the cold ga s in the 
spiral arms to clump together (|Dobbs et al ] l200fj ; iDobbsl 
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Figure 2. The column density of the disc is shown after 240 
Myr for the models with surface densities of 4 (top) 10 (middle) 
and 20 Mq pc~~ 2 (lower). The structure is similar for all surface 
densities, although more high density gas is apparent in the arm 
and interarm regions for the 10 and 20 Mq pc — 2 cases. The scales 
have been chosen so that they are proportional to the average 
column density. 

2008) and these clumps emerge from the spiral arms as dis- 
tinct spurs. 

In Fig. 2, the results for the different surface densities 
are displayed. The corresponding time is 240 Myr, and E = 
4 (top), 10 (middle) and 20 (lower) M Q pc" 2 . Even the 
lower surface density case shows a high degree of substruc- 
ture, since there is still a large component of cold gas. In all 
cases, since the cold gas is the most dense, this dominates 



Figure 3. The pressure is plotted against number density for 

100.000 random particles. The distribution shows gas at stable 
points corresponding to temperatures of 100 (n~20 cm -3 ) and 
10 4 K (n~0.1 cm -3 ) with unstable gas between these regimes. 

the structure. There is little difference between the higher 
surface density cases. However there is comparatively more 
dense gas, both in the spiral arms and the interarm struc- 
tures, than in the 4 Mq pc -2 run. For the higher surface 
density discs, more of the gas reaches very low tempera- 
tures, of ~ 10 K. This very cold gas is denser, and remains 
at high densities for longer. The density is also noticeably 
lower in the interarm regions compared to the lower surface 
density case. 

3.1.1 Thermal distribution of the gas 

In the previous section, we highlighted that cooling of gas in 
the spiral shocks is sufficient for gas to form dense clumps 
in the spiral arms, which in turn lead to interarm features. 
Here we discuss in more detail the thermal distribution of 
the gas. Overall, except for the lowest surface density case, 
the largest proportion of gas ends up in the cold phase. 

Figure 3 shows the pressure distribution versus number 
density of a selection of particles from the disc. Generally 
there is a tight correlation between the pressure and den- 
sity. The distribution corresponds to a two phase medium, 
with stable points at ~ 10 4 K and 100 K. There is however 
gas with 0.01 < n < 1 cm -3 which exhibits a lower than 
expected pressure. This represents gas which has recently 
passed through the spiral arms, but has not yet heated back 
up to equilibrium pressure. This gas is only a few per cent 
of the total amount of gas. Conversely, there are no gas par- 
ticles which have higher than expected pressures, since we 
do not include any heating of the dense gas by supernovae. 

In Fig. 4 we show the distribution of temperature for 
the different surface density simulations. All show peaks at 
around 1 00 and 10 4 K, as expected from the two phase 
model of iFieldl (|1965| ). There is more warm (10 4 K) gas 
in the lower surface density models, whilst at higher sur- 
face densities, more of the gas cools to temperatures be- 
low 100 K. There is also a significant amount of gas in 
the thermally u nstable regime, as has been noted in previ- 
ous simulations (|Gazol et al.ll200ll ; iKritsuk fc Normanll2002l ; 
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Figure 4. The distribution of temperatures for a surface density 
of 4 (dotted), 10 (solid) and 20 (dashed) M pc -2 . This distri- 
bution corresponds to a time of 240 Myr. Peaks are evident at 
100-200 K and 10 4 K, the stable points for cold and warm phases 
of the ISM respectively. 



100 




100 150 
Time (Myr) 



250 



Figure 5. The percentage of gas which is cold (T sC 150 K) 
(blue), warm (T 5000 K) (red) and in the unstable regime (150 
K sC T 5000 K) (green) is shown versus time for the 4 (dotted), 
10 (solid) and 20 (dashed) Mq pc -2 surface density calculations. 



Vazqu ez-Scmadcni ot al. 2003; d e Avillez fc Breitschwerdtl 
2004 iPiontek fc Ostrikerll2005h . In the 10 M pc" 2 model 
plotted here, the photodissociation length scale is l p h = 35 
pc; there are only small differences in the distribution for the 
other length scales, mainly at low temperatures where the 
amount of molecular gas, and therefore the mean molecular 
weight fi differs. 

The evolution of the amount of gas in the cold (T ^ 
150 K), unstable (150 K < T < 5000 K) and warm 
(T ^ 5000 K) phases is shown in Fig. 5. Clearly the amount 
of cold gas in the disc increases with surface density. At the 
outset, all the gas is in the warm regime, at 7000 K. Gas 
which does not immediately enter a spiral arm cools to a 



s 


4 Mq pc~ 2 


10 Mq pc~ 2 


20 Mq pc" 2 


Cold 


40% 


70% 


85% 


Unstable 


46% 


25% 


13% 


Warm 


14% 


5% 


2% 



Table 2. The amount of gas in the cold, unstable and warm 
phases is shown after 250 Myr for the three different surface den- 
sity simulations. 



few 1000 K, in the thermally unstable regime. Gas entering 
the spiral shock cools to much lower temperatures of a few 
100 K. As the shock develops further, the gas continues to 
cool, to temperatures of ~100 K. From Fig. 5, we see that 
the amount of cold gas grows over a period of approximately 

30 Myr, comparable to the typical crossing time in the spiral 
arms. Much of the gas cools to temperatures of < 150 K, 
but a small proportion of the gas also heats back up to the 
warm stable regime. After further spiral arm passages the 
amount of gas in each phase does not change significantly. 

Over the course of each simulation, the evolution of the 
cold and unstable phases appear synchronised, whilst the 
amount of warm gas is fairly constant after 100 Myr. This 
suggests that transitions occur predominantly from the cold 
to unstable phases and vice versa. However, we should again 
stress that we do not include heating from supernovae, or 
from other forms of stellar feedback (e.g. winds) in these 
calculations. 

By 200 Myr (approximately 2 rotations), the amount of 
gas in each phase is relatively stable. Except for the low sur- 
face density case, the largest proportion of the gas is cold. 
Table 2 compares the amount of gas in each phase for the 
different simulations after 250 Myr. There are only few ob- 
servatio nal estimates of the frac tions of gas in the different 
regimes. iHeiles fe Trolandl (|2003l ) estimate that ~ 40% of the 
HI in the solar neighbourhood is cold, 30% is unstable and 
30% is warm. In our simulations, however, the cold compo- 
nent includes molecular gas. Thus discounting the molecular 
gas (which is predominantly cold), the percentages of HI in 
the cold, unstable and warm regimes are approximately 62, 

31 and 6% for the 10 Mq pc -2 calculation. 

Several numerical simulations consider cooling and 
heating in local regions of the ISM, usually applying pe- 
riodic boundary conditions. IPiontek fc Ostrikerl (|2005l ) ob- 
tain similar proportions to the 10 and 20 Mq pc -2 results 
(which represent a comparable average n) for 3D simulations 
of MRI driven turbulence, as do lKim et all (|2008l ) for ID cal- 
culations of gas subj ect to spiral shocks. For 2D simulations 
of driven turbulence, iGazol et all (l200lh find half the gas is 
unsta ble and a quarter warm or cold. lAudit fc Hennebelld 
(|2005l ) also model turbulence and obtain relatively low frac- 
tions of cold and unstable gas (~ 30% total) but this may 
well be due to their boundary conditions; they use inflow 
conditions with a constant injection of warm gas, rather than 
a fixed mass of gas. 

Few global sim ulations have inves t igated the thermody- 
namics of the ISM. IWada fc Normanl (|200lh include stellar 
feedback as well as the thermodynamics in their 2D mod- 
els of galactic discs, although the structure formed in their 
simulations is due to self gravity rather than spiral shocks. 
They find approximately 2/3 of the gas lies in the cold 
regime, 20% is unstable and the remainder is hot gas. Recent 



non-isothermal 3D sim ulati ons of flocculent spirals i nclud e 
iTasker fc Brvanl (|2008h and iRobertson fc Kravtsovl l|2008h . 
but they adopt minimum temperatures of 300 and 100 K 
respectively, and so cannot properly address this question. 

3.1.2 Spatial distribution of warm and cold gas 

Figure 6 displays cross sections of the temperature for a 
section of disc (with z = 0), for the 3 different surface den- 
sities. Evidently there is much more cold gas at higher sur- 
face densities (lower panels), which is situated in the spiral 
arms and the interarm clumps. The warm gas is only visi- 
ble between the cold clumps, although it is fairly ubiquitous 
across the disc. Essentially the warm gas behaves indepen- 
dently of the cold, and from particle plots, the warm phase 
a ppears to shock ear l ier th an the cold gas, as seen in Fig. 2 
of iDobbs fc Bonnelll l|2007h . 

The scale height of the warm gas is much larger than 
the cold (Fig. 7), with warm gas extending to 400 pc or so 
above the plane of the disc. The cold gas is confined within 
100 pc of the plane. It is also apparent from Figs 6 and 7 
that the volume of cold gas (especially below 100 K) is much 
smaller compared to the percentage of cold gas by mass. 

3.2 Evolution of molecular hydrogen 

The remaining part of this paper discusses the formation 
and evolution of molecular hydrogen. The total amount of 
H2 in the disc is compared for the different surface densities 
and assumed photodissociation scale lengths in this section. 
The evolution of molecular hydrogen is also considered in 
Sections 3.2.3 and 3.2.4, in context with the other variables 
in the simulation, i.e. density, temperature, as well as the 
passage of gas through a spiral shock. We present many of 
the results in this section in terms of the H2 fraction, which 
we define as /h 2 = 2n(H 2 )/(n(H) + 2n(H 2 )). 

The time evolution of the percentage of H2 in the disc 
for the different simulations is shown in Fig. 8. As expected, 
the amount of H2 increases with surface density, and the 
scale length for the photodissociation, \ v h which determines 
the H2 column density. From Fig. 8, the percentage of H2 
scales approximately linearly with the surface density. This 
same dependence on surface density was also fo und in pre- 
vious isothermal calculations (|Dobbs et al-l feoOG). From the 
lower panel comparing Z p h we see that a higher estimate 
of N(H2) increases the self-shielding of the molecular gas 
(Eq. 3) and thus allows more H2 to survive. However, the 
dependence of the mean molecular fraction on Z p h is rela- 
tively weak: an increase in Z p h of almost a factor of seven 
increases the mean molecular fraction by no more than a 
factor of two. As realistic values of l p h lie within the range 
15 < Zph < 100 pc, this demonstrates that the uncertainty 
introduced by our simplified treatment of H2 self-shielding 
is unlikely to be very large. 

The fraction of H2 is seen to peak after around 180 Myr, 
corresponding to about 2 spiral arm passages (for gas at the 
midpoint of the disc). The amount of H2 then levels out at 
approximately 7, 20 and 38 % for the 4, 10 and 20 M pc" 2 
surface density discs respectively. For comparison, at the so- 
lar radius, the surface density of H2 from models of the ISM 
is thought to be 1 - 2 Mq pc -2 , whilst the total surface den - 
sity is - 10 M pc" 2 (including He) l|Wolfire et al.ll2003h . 



Cooling in spiral shocks 7 




Figure 6. A cross section (where z = 0) of the temperature is 
displayed for a section of the disc for the 4 (top), 10 (middle) and 
20 (lower) Mq pc" 2 surface density discs. There is more cold gas, 
< 100 K, in the higher surface density calculations. 




x(kpc) 

Figure 7. A cross section showing the temperature in the zx 
plane (y = 0) for a section of the disc. The two areas with cold 
gas (dark blue) on the left and right represent the spiral arms. 
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Figure 8. The percentage of molecular gas in the disc is shown 
versus time. In the top panel, the results for the different surface 
density cases are plotted, whilst the lower panel shows the cases 
for different photodissociation length scales (the assumed distance 
to a B0 star). 



In our E = 10 M© pc - simulation, the surface density of 
H2 is 2 Mq pc~ 2 . Observations suggest the percentage of 
neutral HI which is molecular ov er the whole o f the Milky 
Way is between 20 and 40% (e.g. lFerrierell200ll ). 



3.2.1 Azimuthal and radial distribution of molecular- 
hydrogen 

Figure 9 shows the H2 fraction versus azimuth, for the fidu- 
cial case with a surface density of 10 Mq pc -2 and l p h = 35 
pc. The H2 fraction is calculated by selecting a ring located 
at 7.5 kpc of width 200 pc. The ring is divided into 200 seg- 
ments azimuthally and the H2 fraction is averaged in each. 
At the earlier time of 120 Myr, the H2 is clearly concentrated 
in the spiral arms. This gas has not reached sufficiently high 
densities that molecular gas survives into the interarm re- 
gions. The gas heats up and almost immediately dissociates 
when leaving the arm. However, after 240 Myr, there are 
peaks situated in the interarm regions, in addition to those 
corresponding to the spiral arms. 

In Fig. 10 we show the radial distribution of H2, HI and 
total HI+H2 after 240 Myr. The surface density is calculated 



Figure 9. The molecular gas fraction is shown versus azimuth 
for a ring located at the midpoint of the disc, at times of 120 
Myr (solid) and 240 Myr (dotted). At the earlier time, there is 
cold, dense molecular gas in the spiral arms, but this gas quickly 
warms up and photodissociates on leaving the spiral arms. On 
the other hand, by 240 Myr, gas in the spiral arms has reached 
much higher densities of H2 and large clumps break away into the 
interarm regions, which have a high molecular gas content. 
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Figure 10. The surface density of HI, H2 and total HI+H2 is 
plotted versus radius for the fiducial calculation with £ = 10 Mq 
pc - 2 and l p ^=35 pc. The corresponding time is 240 Myr. The 
surface density profile of HI is relatively flat, whereas that for the 
molecular gas falls off with radius. 



by dividing the mass of gas in a ring at a given radius, by 
the area occupied by the ring. The HI is fairly flat across 
the disc, whilst the H2 falls off with radius. The density of 
the gas tends to be underestimated at the edge of the disc 
though, so there is less molecular gas than expected. The 
decreasing H2 surface density occurs since the spiral arms 
are further apart at larger radii and the molecular gas is 
predominantly in, or formed in, the spiral arms. 



3.2.2 PDFs of the gas density 

The (mass- weighted) PDF is shown in Fig. 11, for both the 
total gas in the disc, and the molecular gas. The cumula- 
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Figure 11. The mass of gas is shown versus density for the fidu- 
cial calculation with S = 10 Mq pc — 2 . The mass for the total gas 
and the molecular gas is plotted, and dashed lines indicate the 
cumulative mass of gas. The mass distribution for the total gas is 
approximated between n=10 cm -3 and 100 cm" 3 by a lognormal 
distribution f(u, <r), where fi = 7.5 and a = 2.46. 

tive mass is also plotted. We find that the total gas PDF is 
approximated by a lognormal distribution, as discussed in 
recen t simulations l)Wada fc Normanll2007l : lTasker fc Brvanl 
2008). However the molecular gas does not follow a lognor- 
mal distribution. The figure indicates that there is negligible 
H2 below 5 cm -3 (the self shielding limit, see Section 3.2.3) 
and thus the total amount of H2 above this density is the 
same as the total H2 in the disc. 



3.2.3 Self-shielding and the formation of molecular 
hydrogen 

Whether gas in the ISM is primarily atomic or molecu- 
lar i s expected to depend on the efficiency of self shielding 
(e.g.lFederman et al.lfl979l : lElmegreenlll993l : iBrowning et all 
2003). We plot the molecular gas fraction versus density in 
Fig. 12. The most striking feature of these results is the large 
increase in molecular gas fraction at a specific density. For 
the fiducial case with Z p h = 35 pc (top panel), this density 
is ~ 5 cm" 3 . This is the density at which self shielding of 
the H2 becomes effective, and the fraction of molecular gas 
can continuously increase with time. Below this density, the 
H2 fraction of the gas is in equilibrium. The time to ob- 
tain molecular fractions of 0.1-1% is very rapid, comparable 
to the timestep of our calculations (0.5 Myr). Hence there 
are few particles with molecular fractions between 10 -6 and 
10 -3 . The colour scheme indicates the mass of gas at a par- 
ticular H2 fraction for a given density. As will be shown 
later, the particles exhibit various trajectories on this plot. 
Gas which increases in density sufficiently to become signifi- 
cantly molecular follows a very similar path from low to high 
H2 fractions, hence the large mass of gas between n = 0.1 
and 100 cm" 3 . 

Figure 13 also compares the H2 fraction for the differ- 
ent simulations. Here the molecular gas fraction is averaged 
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Figure 12. The fraction of molecular gas is plotted against num- 
ber density for all the gas in the disc. The colour scheme indicates 
the amount of mass with a particular H2 fraction for a given den- 
sity. 

over all particles at a given density, although the spread in 
H2 fraction is similar to Fig. 12 for each case. The top panel 
compares the simulations with different surface density. The 
dependence of the H2 fraction on density does not change 
for the different surface density calculations. Clearly when 
E = 20 Mq pc" 2 though, much more of the gas resides at 
high densities, and therefore high H2 fractions. The situ- 
ation changes when considering different Z p h values (lower 
panel) . The effect of increasing Z p h is to decrease the thresh- 
old density at which self shielding becomes important. This 
threshold depends on N{H2) = n(Hi) Z p h rather than n. So 
when Z p h = 100 pc the gas only has to exceed densities of 
n ~ 2 cm -3 to start attaining a high molecular gas fraction 
(and as seen in Fig. 8, lower panel, there is more molecular 
gas overall), compared to n ~ 8 cm -3 when Z p h = 15 pc. 

Figure 13 indicates that a large fraction of the gas is 
molecular at moderate densities, e.g. 10-30% for n — 2 — 10 
cm" 3 and 50% or more once the density exceeds 100 cm -3 . 
We stress again that this is a mean fraction, which includes 
gas that is currently decreasing both its density and H2 frac- 
tion, possibly having already spent many Myr at densities 
above the self shielding limit. Furthermore the H2 fraction 
varies with metallicity, extinction and UV flux, all of which 
are assumed to be constant in our model. For c ompa rison, 
in calculations performed by IPelupessv et al l (|2006h . the 
molecular gas fraction at n = 100 cm -3 is ~ 0.4, 0.9 or 
> 0.9, depending on met allicity, UV flux and the form ation 
rate of H2 on grains. In iKovama fc Inutsukal (|2000l ). sub- 
stantial molecular fractions are not reached until n = 1000 
cm" 3 , but this is probably a consequence of lower column 
de nsities, since they consider a localised slab of gas. Finally, 
m Gl over fc Mac Lowl (|2007bh , turbulent clouds with mean 
densities of 10 and 100 cm -3 develop molecular gas frac- 
tions of ~ 0.1 and 0.4 respectively within only a few Myr, 
consistent with the fractions we find here. 

We can also consider the evolution of particles as their 
density and molecular gas fraction change. Figure 14 shows 
evolutionary tracks of 3 particles from the simulation with 
Z P h = 35 pc and E = 10 Mq pc -2 . The top panel is the most 
typical (i.e. frequent) path for a particle. As a particle enters 
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Figure 13. The molecular gas fraction is shown averaged over 
particles at a given number density. In the top panel, results are 
plotted from the different surface density calculations. The H2 
fraction at different densities is similar regardless of the surface 
density. The lower panels compare results with different photodis- 
sociation approximations. When assuming a large scale length 
and therefore column density, the threshold for self shielding shifts 
to lower densities (100 pc case). 



a spiral shock, its density increases, and the particle evolves 
along the right hand side of the distribution of particles in 
Fig. 12. At some point, the density of the particle starts to 
decrease, and the particle evolves to the left in H2 fraction 
versus n space. This accounts for the large spread in the 
H2 fraction apparent for densities > 5 cm" 3 in Fig. 12. But 
although the particle's density decreases, it still retains a 
high molecular gas fraction. Thus the evolution in H2 frac- 
tion is quite different depending on whether the particle's 
density is increasing or decreasing. The reason for this is 
because once the fluid element has attained a high molec- 
ular gas fraction, the column density N(Hz) is large and 
self shielding (denoted by /shield in Eq. 2) prevents the H2 
from rapidly becoming photodissociated. Thus the gas can 
exhibit higher H2 fractions at lower densities compared to 
before the gas enters the shock (where N(H2) and the degree 
of self shielding is negligib le) . A similar scenario is described 
in iPelupessv et al.l |2006) , who even find that molecular gas 
can survive the transition from the CNM to the WNM. How- 



ever in reality many of the trajectories in Fig. 12 would be 
interrupted at high densities by star formation processes. 

The lower panels in Fig. 14 indicate the evolution of two 
gas particles which behave somewhat differently. The mid- 
dle panel indicates an example where a particle exhibits a 
high molecular gas fraction for most of the calculation. This 
is because once the gas attains a large molecular fraction, it 
encounters another spiral shock before becoming photodis- 
sociated to a very low H2 fraction. The gas particle in the 
lowest panel exhibits the opposite behaviour. It lies in the 
warm regime, at a larger scale height. The gas does not ex- 
perience a shock since it is too warm and diffuse. Thus the 
density does not become sufficiently high for the fluid ele- 
ment to attain a significant H2 fraction. 



3.2.4 Evolution of physical variables during the simulation 

This section describes how the density, H2 fraction and tem- 
perature evolve with time for specific gas particles in the 
simulation. Fig. 15 shows the evolution of the various physi- 
cal quantities with time, as well as the distance to the near- 
est potential minimum, for the three particles already de- 
picted in Fig. 14. The corresponding galactic radii for the 
particles are 7.4 kpc (top), 5.1 kpc (middle) and 6 kpc (lower 
panel). For particles at larger galactic radii, there are corre- 
spondingly fewer minima in the variation of the distance to 
the potential minimum with time, indicating fewer passages 
through the stellar spiral arm. 

The top panel shows the typical behaviour of gas as it 
travels round the disc. After 30 to 40 Myr, the gas parti- 
cle passes the potential minimum and enters a spiral shock. 
The gas experiences a sharp increase in density, and a cor- 
responding drop in temperature to ~ 100 K. The increase 
in molecular gas fraction is similarly rapid, changing from 
~ 10" 6 to ~ 10~ 3 in a single timestep (0.5 Myr). A further 
increase to the peak fraction of ~ 0.1 takes around 10 Myr. 
The gas then remains in the spiral arm for another 40 to 
50 Myr. After leaving the spiral arm, the gas becomes less 
dense and heats up. At the same time, the gas experiences a 
sudden drop in the H2 fraction, decreasing to ~ 10 -6 in 10 
Myr or so. The gas enters a spiral shock for the second time 
after 190 Myr. By this time, the spiral shock is more devel- 
oped and the gas reaches higher densities. Consequently it 
cools to ~20 K and attains a molecular gas fraction of 0.7. 
Again the time to acquire this fraction is short, around 10 
Myr. At this point in the simulations there is much more in- 
terarm structure, with clumps breaking off into the interarm 
regions. The gas retains a higher density and H2 fraction for 
a longer period, including time when gas has left the spiral 
shock, but drops down to lower densities and a very small 
H 2 fraction by 300 Myr. 

The middle panel in Fig. 14 shows the evolution of gas 
which retains a high density throughout most of the simu- 
lation (as also described in Section 3.2.3). This gas is nearer 
the inner edge of the disc and experiences more frequent 
spiral arm passages. The density does not decrease suffi- 
ciently for the H2 fraction to return to 10~ 6 . Consequently 
the molecular gas fraction stays between 10~ 3 and 0.1, whilst 
the temperature remains around 100 K (thus the gas remains 
in the cold regime for the rest of the simulation). Finally in 
the third panel, the gas never attains a significant fraction 
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Figure 14. The trajectories of three particles are shown in H2 
fraction versus density space. The colour of the plotted lines 
changes with time, with red at t = Myr to black at 310 Myr 
(arrows also indicate the direction the particles evolve over time) . 
In the top panel, the gas (twice) reaches a high molecular gas 
fraction and then drops down to lower densities with very low 
molecular gas fractions. This pattern is typical for most particles 
in the disc. The average galactic radius for each of the particles is 
also indicated on the panels. The evolution of H2 fraction, density 
and temperature are shown with time in the corresponding plots 
in Fig. 15. 



Figure 15. The evolution of density, H2 fraction, temperature 
and column density is shown versus time for 3 particles (a different 
particle in the top, middle and lower panels). The distance to the 
nearest potential minimum was also calculated and is shown in 
these figures. The top panel shows fairly typical behaviour as 
the gas passes through the spiral shock, where the H2 fraction 
increases in the spiral arm and then decreases. In the middle 
panel, after the first spiral shock, the gas retains a high density 
and H2 fraction throughout the simulation. For the third panel, 
the gas density never gets very large, the gas has negligible H2, 
and remains in the warm phase of the ISM. The galactic radii of 
the particles are indicated on each of the panels. 
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of H2 . The gas remains too warm and diffuse throughout the 
simulation. 

We compared the profiles of 100 random particles, find- 
ing that the majority of particles (60-70%) exhibit evolu- 
tionary patterns similar to the top panel of Fig. 14, i.e. the 
molecular gas fraction increases and decreases by orders of 
magnitude during and after each spiral passage. A few parti- 
cles have high densities and molecular gas fractions for most 
of the simulation, similar to the middle panel, but these tend 
to be nearer the centre of the disc. A particular gas particle 
does not always exhibit the same behaviour through each 
spiral arm, e.g. gas may be low density and not molecular 
at the beginning of the simulation, but increase its peak 
density during successive spiral arm passages. Conversely a 
gas particle may be dense and have a high molecular gas 
fraction before 100 Myr, but then move into a low density 
region between spiral arms, and become low density atomic 
gas for the rest of the simulation. 

3.2.5 Timescales for the formation and destruction of H2 

The questions of how quickly molecular clouds form, and 
how long they live, are central to any model of star forma- 
tion. The simulations performed for this study can provide 
an insight into these timescales. 

Figure 16, shows the time taken for the gas to achieve 
an H2 mass fraction of 0.5. The characteristic time required 
for raising the molecular fraction from 0.1 to 0.5 is around 
10 Myr, which is roughly the time needed to cross the spiral 
arm. We stress here that the values given in Fig. 16 are up- 
per limits, since we do not have the resolution to treat the 
intern al turbulent motions in the spiral arms fully consis- 
tently. iGlover fc Mac Low! (|2007bh demonstrated that tur- 
bulence at and below the GMC sc ale can induce H2 forma - 
tion on an internal crossing time. iDobbs fc Bonnelll (|2007T l 
showed that when structured gas is passed through a shock 
it generates turbulence that is consistent with t he typical ve- 
locity scaling laws found in molecular clouds (|Larso n 1981; 
iHever et al.ll2004r ). The fact that the typical timescale for 
going between molecular fractions from 0.01 to 0.1 is faster 
than from fractions 0.1 to 0.5 is indicative of the resolu- 
tion constraint our calculations. It is therefore likely that 
the unresolved velocity structure in the spiral arms would 
accelerate the formation of H2 to rates higher than we quote 
here. Self-gravity would increase these rates further still and 
as such, one would expect the peaks in the formation time 
distributions to move to lower values. 

The ability of a parcel of gas (in this case, an SPH 
particle) to gain a large H2 fraction, depends on its path 
through the spiral arm, and hence the level of compression 
that it experiences. Two effects are important here. The first 
is 'orbit crowding', whereby gas streams from different radii 
are brought together by the spiral potential due to their 
relative motions. This increases the density of the gas that 
undergoes a shock during its passage through the spiral. The 
second effect is simply that gas entering the spiral arm will 
undergo a stronger shock if the gas already residing in that 
part of the spiral arm is cold. 

At earlier stages in the simulation, the spiral shock is 
not well developed. Warm gas experiences compression and 
achieves low (~100 K) temperatures, but usually only mod- 
erate H2 fractions (e.g. Fig. 9, Fig. 15 top panel). However, 



subsequent gas entering the shock now encounters cold gas. 
This gas experiences not only compression by the shock, 
but collides with cold clumps in the shock. As a result, the 
gas experiences stronger compression, cools to lower tem- 
peratures and can obtain a higher H2 fraction. With time, 
this process repeats, and the gas in the spiral shock be- 
comes colder and more molecular. Eventually a steady state 
is reached (Fig. 8), by which time gas can enter the shock 
and achieve high H2 fractions in 10 to 20 Myr. The formation 
of GMCs is thereby facilitated by the pre-existence of d ense, 
cold material, as hypothesised bv lPringle et al.l (|200ll ). 

The density reached in the spiral arms also depends on 
the density and temperature of gas entering the shock (e.g 
Fig. 15, middle and lower panels). Gas which is in the mid- 
plane of the galaxy, or which has not fully heated up after 
passing the previous arm, is denser and colder and is more 
likely to achieve a high H2 fraction (Fig. 15, middle). As 
discussed in Section 3.3.1, much of the gas lingers at densi- 
ties around 0.1 to 1 cm -3 and temperatures of several hun- 
dred Kelvin. This places it in the regime where self-shielding 
starts to become important (see Figures 11 & 12). As such, 
the gas remains in a regime where it can easily increase its 
molecular fraction during the compression produced by the 
spiral arm passage. Gas at low density, e.g. above the mid 
plane, may pass through the spiral arm without obtaining a 
significant increase in H2 fraction (Fig. 15, lower). Overall, 
the combination of orbit crowding and collisions between 
parcels of cold gas that is described above creates large cold 
molecular complexes. Between passages through the spiral 
arms, these parcels of cold gas typically do not maintain a 
high H2 fraction. 

Finally Fig. 17 indicates the time gas spends above 
given molecular fractions. The gas tends to spend more 
time at higher H2 fractions. We can relate this behaviour 
to Fig. 14, which shows the trajectories of particles in H2 
versus n space. The particles' trajectories form loops, which 
increase in size as the gas reaches higher densities. Gas which 
reaches the high density regime is highly molecular and takes 
a longer time to photodissociate back to lower H2 fractions. 
Fig. 17 indicates that gas which reaches an H2 fraction of 
0.01 usually retains this fraction for only 10 to 20 Myr. How- 
ever gas which goes on to obtain H2 fractions of 0.5 retains 
this H2 fraction for ~40 Myr. This is clearly an upper limit 
for the expected cloud lifetime, as we do not include star 
formation and stellar feedback. The strong UV radiation 
and winds from newly formed massive stars would acceler- 
ate molecular cloud disruption. 



4 DISCUSSION 

The calculations presented in this paper demonstrate 
that the atomic and m olecular cooling discussed in 
Gl over fc Mac Low! l|2007al ) permits the formation of molec- 
ular clouds in spiral shocks. Warm gas flowing into the spiral 
potential is compressed, enabling it to cool to low temper- 
atures, and initiating the rapid formation of H2 in the cold 
gas. By 'rapid' we mean that large quantities of H2 form in 
roughly a crossing time, consistent with the ideas presented 
by a number of previous studies (jVazauez- Scmadcn i et all 
1996 ; [Ballesteros-Paredes et al.l 1 19991: iKovama fc Inutsukal 



200ol : lElmegreenll2000l : iPringle et alj|200ll) 
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Figure 16. This histogram gives the distribution of timescales 
over which the gas reaches certain molecular gas fractions. The 
timescales denote the time for the H2 fraction of a particle to 
increase from 0.001 to 0.01, 0.01 to 0.1 and 0.1 to 0.5, as indicated. 
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Figure 17. This figure shows the distribution of timescales gas 
remains over molecular gas fractions of 0.01, 0.1 and 0.5. Gas 
typically spends about 40 Myr with a molecular gas fraction of > 
0.5. This is however an upper limit on a molecular cloud lifetime, 
since stellar feedback is not included in these models. 



Our results differ from the two competing theories 
of GMC formati on. In one, cloud formation is driven by 
self-gravity (e.g. lElmegreenl Il979). and H2 formation oc- 
curs on the cloud collapse timescale (which may be signifi- 
cantly longer than the free-fall time for clouds that have a 
high level of turbulent or magnetic support). In the other, 
GMC formation occurs through the coalescence of many 
small cloudlets, wh ich are dense but not self-gravitating 
l|Pringle et al.| [2001). Our calculations cannot directly ad- 
dress the first of these scenarios, since we do not include the 
effects of the self-gravity of the gas. However, our results 
do highlight the importance of the gravitational potential 
of the spiral arm for driving the formation of GMCs. Orbit 
crowding within the spiral arm leads to cloud coalescence, 



as discussed in Section 3.2. 5. The final p i cture is therefore 
similar to that proposed bv lPringle et al.l l|200ll ). 

The rapid formation of molecular clouds as 
suggested in the current study, as well as in a 
number of more det ailed chemical and d y nami- 
Ballcstcro s-Paredes et al. Il999l; 
20071; iGlover fc Mac Low! 



cal 



more 
calculations (e.| 



Heitsch ct al 



2006 



2008 



Ballcstcros-Parcdcs & Hartmann 
| 2007allbl; lHartmann et all |200ll; 

Hennebelle et al.l 120071 . 120081 ; IVazquez-Semadeni et"aH 
20061 . 120071 ). agrees well with current observational data. 
For example, it is consistent with the geometrically derived 
fast star formation timescales in 14 nearby spiral galaxies 
observed with the VL T and the Spitzer infrared satellite. 
iTamburro et all (|200&t ) measure the average angular offset 
between the Hi and 24 /jm emissivity peaks along spiral 
arms and find typical timescales for the onset of star 
formation in this sample of below 5 million years. Earlier 
studies on the position of dust emission, Hi, Hn, and Ha 
line from young stars measured across spiral shocks in 
M51 and M83 dAllen et al.l Il98rj ; iTilanus fc Alien] Il989l ; 
iLeisawita et al] ll989T l" come to similar conclusions. Short 
cloud formation timescales are also consistent with the 
stellar age spread in molecular clouds which exhibit a short 
peak of star-formation activity of 1 to 2 mill ions years 
duration l)Hartmann et al1l200ll : lHartm ann 2003), possibly 
preceded by an acceleration period with very low a ctivit y 
of up to 10 million years (e.g. lPalla fc Stahlerllioool . I2002T I. 
Altogether, we conclude that a cloud formation timescale 
of about 10 million years or less as suggested by our 
calculations is consist ent with observational data (see also 
lElmegreenl [2OO0I . I2007T I . 

Our study addresses another fundamental point of 
GMC formation: local self-gravity is not required to form 
molecular clouds. The 'clouds', or complexes of molecular 
gas, are simply a result of a change in chemical state brought 
about by the processes described above. It is not necessary to 
assume that GMCs are dominated by self-gravity. In partic- 
ular, GM Cs need not be in v irial equilibirum, as previously 
noted bv lPringle et~aD (|200lh . The path to molecular cloud 
formation is already underway at densities of 0.1 to 1 cm -3 
and at temperatures of around 100 - 300 K, the regime in 
which the self-shielding becomes effective. The free-fall time 
at this density is around 50 Myr. Since this is longer than the 
typical cloud formation timescale, it is unlikely that grav- 
ity can play a dominant role in the formation of molecular 
clouds. Additionally, self-gravity has to overcome turbulence 
present in the ISM, and tidal shear from the galaxy poten- 
tial. Given these conditions, it seems likely that globally 
unbound GMCs can form, in which the star formation effi- 
ciency would be naturally quite low l|Clark et al.ll2005h . 



5 CONCLUSIONS 

The simulations presented in this paper indicate that cooling 
in spiral shocks leads to the formation of molecular clouds 
in the spiral arms of a galaxy. Self gravity of the gas is not 
required for this process. Cooling of the gas is very quick as 
gas enters the shock (almost instantaneous), hence cold gas 
clumps can agglomerate into larger structures, as described 
in previous i sothermal simulations l)Dobbs fc Bonnellll2006l : 
iDobbs et ail 120061 ; iDobbsl lioOsT ). within the spiral shock. 
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Larger GMCs form after successive shocks, once both the 
arm and interarm material has become colder and denser. 
Again, similar to previous calculations, these clumps are 
sheared into interarm spurs. 

For a disc of surface density comparable to the so- 
lar neighbourhood, around 70% of the gas lies in the cold 
regime, and 20% is molecular. The transition from atomic to 
molecular gas is shown to be strongly dependent on the self 
shielding limit. Spiral shocks compress the gas sufficiently 
that large amounts of gas exceed the self shielding limit, 
thus enabling rapid formation of H2. Gas obtains moderate 
H2 fractions immediately after compression. The maximum 
H2 fraction obtained then depends on the path of the gas 
through the shock. Gas which encounters little high density 
material may obtain a fraction of 0.01 or so for a short period 
(10-20 Myr) after which photodissociation reduces the H2 
fraction. On the other hand gas may encounter dense, cold 
material, in which case it will undergo a stronger shock, be 
subject to greater cooling and become predominantly molec- 
ular. The time for the transition to fractions of 0.5 or so 
is typically 10 to 20 Myr, although this is a strict upper 
limit because our calculations cannot fully resolve internal 
turbulence and neglect self-gravity of the gas. True cloud 
formation timescales therefore will be considerably shorter. 
Without stellar feedback however, this gas tends to remain 
at high densities, and therefore retain appreciable H2 frac- 
tions for substantially longer times (~ 40 Myr). We note the 
cloud formation time is generally not a well defined quan- 
tity, and suggest that determining the time from which gas 
enters a spiral arm, to the point at which it becomes ob- 
servable (in CO), may be a more robust measure in future 
calculations. 



from http://www.astro.ex.ac.uk/people/dprice/splash/ 

i|Pricdl2007n ~ 



APPENDIX A: RESOLUTION STUDY 

To investigate the effects of resolution, we repeated our fidu- 
cial calculation, with E= 10 M© pc -2 and Z p h = 35 pc using 
1 and 4 million particles. The results for the 4 and 8 million 
particle calculations are very similar, although some discrep- 
ancies arise, particularly regarding the amount of cold and 
molecular gas, when there are only 1 million particles. 

On global scales, there is barely any difference in the 
structure of the disc for all the different resolution simula- 
tions. To highlight possible differences, Fig. Al zooms in on 
a section of the disc. The structure is very similar for the 4 
and 8 million particle calculations. For the 1 million parti- 
cle simulation, there is slightly less structure both along the 
spiral arms, and in the interarm spurs, since the small scale 
structure is not as well resolved. 

We also compared the thermal distribution of the gas, 
and the amount of H2 formed at different resolutions. Fig- 
ure A2 indicates that the thermal distribution is very sim- 
ilar at different resolutions. The distributions mainly differ 
at low temperatures where the cold, dense regions are not as 
well resolved at lower resolutions. We compared the amount 
of molecular gas after 240 Myr in the 1, 4 and 8 million 
particle simulations, finding that 12, 21 and 22 % of the gas 
in the disc is molecular. Thus with only 1 million particles 
the amount of molecular gas is significantly underestimated. 
Again the amount of molecular gas is similar in the 4 and 8 
million simulations, although the 22% for 8 million particles 
still represents a lower limit. 
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